Asymptotic approximation for the solution to a semilinear parabolic problem in a thin star‐shaped junction

A semilinear parabolic problem is considered in a thin 3‐D star‐shaped junction that consists of several thin curvilinear cylinders that are joined through a domain (node) of diameter The purpose is to study the asymptotic behavior of the solution u_ε as ε→0, ie, when the star‐shaped junction is transformed in a graph. In addition, the passage to the limit is accompanied by special intensity factors and in nonlinear perturbed Robin boundary conditions. We establish qualitatively different cases in the asymptotic behavior of the solution depending on the value of the parameters {α_i}and {β_i}. Using the multiscale analysis, the asymptotic approximation for the solution is constructed and justified as the parameter ε→0. Namely, in each case, we derive the limit problem (ε=0)on the graph with the corresponding Kirchhoff transmission conditions (untypical in some cases) at the vertex, define other terms of the asymptotic approximation and prove appropriate asymptotic estimates that justify these coupling conditions at the vertex, and show the impact of the local geometric heterogeneity of the node and physical processes in the node on some properties of the solution.     

Asymptotic properties of a stochastic chemostat including species death rate

This paper deals with a stochastic system which models the population dynamics of a chemostat including species death rate. On the basis of the theory on Markov semigroup, we demonstrate that the probability densities of the distributions for the solutions are absolutely continuous. The densities will convergence in L¹ to an invariant density or weakly convergence to a singular measure under appropriate conditions. We also give the sufficient criteria for extinction exponentially of the species. To be specific, when D₁>D and the strength of perturbation is relatively small, we derive a precise threshold for the species survival.     

One‐level and two‐level finite element methods for the Boussinesq problem

Previous works on the convergence of numerical methods for the Boussinesq problem were conducted, while the optimal L²‐norm error estimates for the velocity and temperature are still lacked. In this paper, the backward Euler scheme is used to discrete the time terms, standard Galerkin finite element method is adopted to approximate the variables. The MINI element is used to approximate the velocity and pressure, the temperature field is simulated by the linear polynomial. Under some restriction on the time step, we firstly present the optimal L² error estimates of approximate solutions. Secondly, two‐level method based on Stokes iteration for the Boussinesq problem is developed and the corresponding convergence results are presented. By this method, the original problem is decoupled into two small linear subproblems. Compared with the standard Galerkin method, the two‐level method not only keeps good accuracy but also saves a lot of computational cost. Finally, some numerical examples are provided to support the established theoretical analysis.     

Qualitative theory for Volterra difference equations

We investigate l^p boundedness, the topological structure of solutions set and the asymptotic periodicity of Volterra functional difference equations. The theoretical results are complemented with a set of applications.     

B‐spline spectral method for constrained fractional optimal control problems

In this paper, we present a direct B‐spline spectral collocation method to approximate the solutions of fractional optimal control problems with inequality constraints. We use the location of the maximum of B‐spline functions as collocation points, which leads to sparse and nonsingular matrix B whose entries are the values of B‐spline functions at the collocation points. In this method, both the control and Caputo fractional derivative of the state are approximated by B‐spline functions. The fractional integral of these functions is computed by the Cox‐de Boor recursion formula. The convergence of the method is investigated. Several numerical examples are considered to indicate the efficiency of the method.     

Analysis of dynamics in a general intraguild predation model with intraspecific competition

This paper is devoted to studying the dynamical properties of a general intraguild predation model with intraspecific competition. We first investigate the stability of all possible equilibria in relation to the ecological parameters, and then study the long time behavior of the solution. Moreover, we provide a detailed analysis of dynamics of a IGP model with linear functional response and intraspecific competition. Our results show that the impact of the intraspecific competition essentially increases the dynamical complexity of the system.     

计算最优控制辛数值方法

针对最优控制问题（OCP）的辛数值方法研究及应用进行综述。主要涉及内容包括,动力学系统为常微分方程描述的一般无约束、含不等式约束和状态时滞的最优控制问题,微分代数方程描述的一般无约束、含不等式约束和含切换系统的最优控制问题,以及闭环最优控制问题。从间接法和直接法两个求解框架出发,重点介绍本课题组在保辛算法方面的研究工作。在间接法框架下,首先基于生成函数和变分原理,将OCP保辛离散为非线性方程组,再数值求解方程组。在直接法框架下,将OCP保辛离散为有限维的非线性规划问题（NLP）,再数值求解。针对闭环最优控制问题,提出了保辛模型预测控制、滚动时域估计和瞬时最优控制算法。研究表明,保辛算法具有高精度和高效率的特点,在航空航天和机器人等领域有着广泛应用前景和价值。     

Optimal variance stopping with linear diffusions

We study the optimal stopping problem of maximizing the variance of an unkilled linear diffusion. Especially, we demonstrate how the problem can be solved as a convex two-player zero-sum game, and reveal quite surprising application of game theory by doing so. Our main result shows that an optimal solution can, in a general case, be found among stopping times that are mixtures of two hitting times. This and other revealed phenomena together with suggested solution methods could be helpful when facing more complex non-linear optimal stopping problems. The results are illustrated by a few examples.     

基于新巴塞尔协议监管下保险人的均值-方差最优投资-再保险问题

本文研究了均值-方差优化准则下,保险人的最优投资和最优再保险问题.我们用一个复合泊松过程模型来拟合保险人的风险过程,保险人可以投资无风险资产和价格服从跳跃-扩散过程的风险资产.此外保险人还可以购买新的业务(如再保险).本文的限制条件为投资和再保险策略均非负,即不允许卖空风险资产,且再保险的比例系数非负.除此之外,本文还引入了新巴塞尔协议对风险资产进行监管,使用随机二次线性(linear-quadratic,LQ)控制理论推导出最优值和最优策略.对应的哈密顿-雅克比-贝尔曼(Hamilton-Jacobi-Bellman,HJB)方程不再有古典解.在粘性解的框架下,我们给出了新的验证定理,并得到有效策略(最优投资策略和最优再保险策略)的显式解和有效前沿.